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Abstract: We describe a few properties of the XXX spin chain with long range interac- 
tion. The plan of these notes is : 
| 1 — The Hamiltonian. 

2 — Symmetry of the model. 

3 — The irreducible multiplets. 
Q-i 

CD ! 4 — The spectrum. 

5 — Wave functions and statistics. 

6 — The spinon description. 
^ | 7 — The thermodynamics. 



Introduction. The XXX spin chain with long range interaction is a variant of the 
spin half Heisenberg chain, with exchange inversely proportional to the square distance 
between the spins. It possesses the remarkable properties that its spectrum is addi- 
tive and that the elementary excitations are spin half objects obeying a half-fractional 
statistics intermediate between bosons and fermions. In this sense, it gives a model 
for an ideal gas of particles with fractional statistics. The model is gapless; its low 
energy properties belong to the same universality class as the Heisenberg model, and 
are described by the level one su(2) WZW conformal field theory. 
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1 The Hamiltonian. 



The Hamiltonian of the trigonometric isotropic spin chain with long range interaction 
is given by g, || : 

where Pij is the operator which exchange the spins at the sites i and j. We restrict 
ourselves to the su(2) case, in which case the spin variables can only take two values: 
Oi = ±. The sum is over all the distinct pairs of sites labeled by integers i, j, ■ ■ ■ ranging 
from 1 to N. 

The spectrum of ([[]), which has been conjectured by Haldane ||, possesses a remark- 
able additivity property as well as a rich degeneracy. It can be described as follow. To 
each eigenstate multiplet is associated a set of rapidities {m p } which are non- consecutive 
integers ranging from 1 to (N — 1). The energy of an eigenstate with rapidities 

{m p } is: 

H\{m p }) = (X>( m p)) \i m p}) with < m ) = (jj) m(m-N) (2) 

The degeneracy of the multiplet with rapidities {m p } is described by its su(2) represen- 
tation content as follows. Encode the rapidities in a sequence of (N — 1) labels or 1 
in which the l's indicate the position of the rapidities; add two O's at both extremities 
of the sequence which now has length (N + 1). Since the rapidities are never equal 
nor differ by a unit, two labels 1 cannot be adjacent. A sequence can be decomposed 
into the product of elementary motifs. A motif is a series of Q consecutives O's, and 
it corresponds to a spin representation of su{2). The representation content of a 
sequence is then the tensor product of its motifs. 

The degeneracy of the spectrum can also be described in terms of path, in a way 
surprisingly similar to the path description of the six- vertex corner transfer matrix B. 



2 Symmetry of the model. 

The symmetry algebra responsible for the degeneracy of the model was identified as the 
su(2) Yangian ||. A Yangian is an infinite dimensional associative algebra generated by 
elements T® b , with n a positive integer, and a, b = ± in the su{2) case. These generators 
satisfy quadratic relations which can be arranged into an Yang-Baxter equation by intro- 
ducing the transfer matrix T(x), with matrix elements, T ab (x) = 5 ab + J2n>o x~ n ~ x T^ . 
The commutation relations then take the following form [|], [F| : 

R(x - y) (1 <g> T{x)) (T(y) ® 1) = (T(y) ® 1) (1 ® T(x)) R(x - y) (3) 

The matrix R(x) is the solution of the Yang-Baxter equation given by: R(x) = x + P, 
where P is the permutation operator which exchanges the two auxiliary spaces. The 
transfer matrix was constructed in || . Its expression is : 

N i 

T ab (x) =5 ab +Y Xf( ) (4) 



with Lij = (1 — 5ij)9ijPij, Oij = Zi/zij with = Zi — Zj, and Xf h is the canonical 
matrix \a)(b\ acting on the i th spin only. The transfer matrix (^) form a representation 
of the exchange algebra (|3]) for any values of the complex numbers Zj . The center of the 
su(2) Yangian algebra (^) is generated by the so-called quantum determinant Det g T(x) 
defined by Q: 

Det q T(x) = T__(x - l)T ++ (x) - T_ + (x - l)T+_(x) (5) 

In the representation (§) , the quantum determinant is a pure number for any values of 
the z/s given by : 

with Aat(x) the characteristic polynomial of the N x N matrix with entries 6^: 
A N (x) = det(x - 6). 

The trigonometric spin chain corresponds to Zj = uj^ with oj a primitive iV root of 
the unity For these values of Zj, the tranfer matrix (§) commutes with the Hamiltonian 
([!]). For Zj = lo\ the matrix G can be diagonalized giving the following expression for 
A N : 

SL N + i 

A N (x) = n(s + - j) (7) 

3=1 / 

3 The irreducible multiplets. 

Solving the model consists in finding all the irreducible components of the Yangian 
symmetry algebra and computing the energy in each of these blocks. For the values of 
the z/s induced by the spin chain, Zj = the representation (|J) is reducible. It is 
completely reducible since the transfer matrix is hermitic: t^ 3 < = t b ^. Each irreducible 
sub- representation possesses a unique highest weight (h.w.) vector |A) which is annihi- 
lated by T + _(x) and which is an eigenvector of the diagonal components T±±(x) of the 
transfer matrix : 



r W |A) = ( ( ++W ° )|A) 



Here, t±±(x) are rational functions in x, but not operators. Since the quantum deter- 
minant (|J) take the same value in any of the irreducible block, these two functions are 
related by : 

M£ + l) =t __(,-i) f++w 

A N (x) 

Hence, only one of them, say t (x), is independent. It uniquely characterizes the su(2) 

Yangian representation. We therefore have to compute all the functions t (x) arising 

from the decomposition of the Yangian representation induced by the spin chain, but 



Obviously, the ferromagnetic vacuum |0) = | + +••• + +) is a h.w. vector: the 

corresponding t (x) is one, and the energy is zero. The h.w. vectors in the one- 

magnon sector are |m) = a; mj crj"|ll), with, 1 < m < (N — 1): the corresponding 
eigenvalue is t (x) = Pl p*^ , with P\{x) = (x + ^±l—m), and the one-magnon energy 

is e(m) = m(m — N). 

In order to determine all the highest weight vectors, we decompose the hilbert space 
into subspaces of fixed magnon number. A M-magnon state has M spin reversed : 

l*> = E i^^-^l") (8) 

m,— ,um 

where cr* denote the Pauli matrices acting on the spin located on the site n. By 
construction, the coefficients if) ni ,---,n M of the M-magnon wave functions are symmetric 
in their indices. The wave function coefficients are unspecified for two coincident indices 
ip... t n,-,n,— By convention, we choose these coefficients to be zero. 

Since these indices range from 1 to N, to any M-magnon state is associated a 
symmetric polynomial *&(zi, • ■ ■ , Zm) in M variables of degree less than (N — 1) such 
that ty(u ni , ■ ■ ■ ,u UM ) = ipni,--,n M - In the following, we restrict ourselves to the class of 
magnon states deriving from polynomials of the following form : 

$(zi, ■ ■ ■ , z M ) = Y[(z p - z q ) 2 R(zi, • • • , z M ) (9) 

with R(zi, • • • , zm) a symmetric polynomial of degree less than (N — 2M + 1). This 
class of states does not include all the states of the spin chain but, as we will see, all 
the highest weight vectors are in this class. 

As explained in the Appendix, the operators T (x) and T + _(x) act on this class 

of states. Therefore, the action of these operators on these magnon states induces an 
action on the polynomials. As shown in the Appendix, we find : 

(•<•) = | 1 + E - . Nil r, I Vl/ « "-) «1 ( J) 



p=l x i 2 U P, 



and 




T + _(x)*(z.)= 1 + E ^N+i n *(*i = M.) (11) 
\ ^x+^-DpJ 

Here we have introduced differential operators D p which have recently been proved 



useful in the Calogero-Sutherland model [10, 11, 12, pi. For the following, we also need 



another set of differential operators D p . Both are defined by : 

Dp Z p Zp -\- ^ QpqKpq 

Dp Z p Zp ~\- ^ ' OpqKpq ^ ] GqpKpq (12) 

q>p q<p 

Here, 6 pq = z p /z pq and, the operator K pq exchanges the positions: K pq z p = z q K pq . 



whereas the D p s are not. On the other hand, the differentials D p commute, 



D p ,D q 



0, but the differentials D p s does not: 



D p ,D q 



(Dp — D q )K P q. The sum of the m 



th 



powers of both differentials form two sets of commuting operators. However, they are 
not independent thanks to the following relation: 

p^x-Dp \ x-Dj \ x-D M ) 

In Eq. ([THp it is understood that the operators are acting on functions symmetric in 
their arguments. 

From eq. (|ll|), we learn that the highest weight vectors correspond to polynomials 
^*(z) with R(z) given by : 

M 

R(z x , ••-,zm) = (YI Zp)<l>(zi, zm) (14) 
p=i 

with (f>(zi, • • • , zm) a symmetric polynomial of degree less than (N — 2M). Thus, in the 
M-magnon sector, there are j^^^liy. independent highest weight vectors. 

From eq. (0), we learn that the eigenfunctions of T (x) are the eigenvectors of 

the commuting hamiltonians of the Calogero-Sutherland model, or equivalently of the 
commuting operators D p . The symmetric eigenfunctions ^^ m ^(z.) with, 

£(Sp) n ^ {mp \z) = E< * {mp} (*-)> 
p p 

are polynomials with degree between 1 and (N — 1) if 1 < m p < (N — 1). Hence, 
the M-magnon highest weight vectors |{^ P }), with wave function given by ^/^ mp \z.), 
are labeled by sets of M integers {mi, ■ • ■ , tum}- Due to the Vandermond prefactor in 
eq.(|9|), these integers never coincide nor differ by a unit. Using the factorisation relation 
(|T3"D, one find that the value of t (x) for these highest weight vectors are : 

, , Pi(x + 1) , , . *L, N + l 
t—(x)= 1K J with P 1 ( X ) = Y[( X + ^— -m p ) (15) 

p= i * 

The dimension of the irreducible multiplets are encoded in the transfer matrix eigen- 
values. The eigenvalues t (x) are given by eq. (0). The remaining eigenvalues t ++ (x) 

are computed from the relation (|J). They can also be written in a product form. The 
result is : 



T(x)\{m p })- FlU 



/ P (x+1) 
1 Po(x) 



\{m p }) (16) 



V 



Pi(x) 

with Pi(x) given in eq. (p~5|) , and Po(x) and Pi(x) factorize An(x) : 

A N (x) = P (x) P 1 (x)P 1 (x - 1). (17) 



One can check that the factorization equation (|17]) admits solutions only if the roots of 



Let us decompose the sequence of rapidities {m p } in elementary motifs as explained 
in Section 1. To each motif of length Q, we associate a canonical transfer matrix defined 
by : 

cab 

T£ oUf (x) = 5 ab + (18) 

1 X — Xq 

where S ab are the matrices forming the spin representation of sit (2) and x is 
the position of the most left label of the motif. It is easy to check that the matrix 
(|18"D satisfy the commutation relations (|3]) . The representation induced by the transfer 
matrix (|16D is then seen to be equivalent to the irreducible tensor product of the transfer 
matrices associated to each motifs: 

N +1 

T(x)= (g) T rnotlf (x + ^— ) (19) 

motifs ~ 

This is proved by comparing the eigenvalues of the diagonal transfer matrix elements 
on the h.w. vector. Therefore, we find that the multiplet of a rapidity sequence is the 
tensor product of each of its motifs. 

We have found one (and only one) irreducible representation for each rapidity se- 
quence. Their direct sum is a vector space of dimension 2^. Therefore, it fills the 
Hilbert space of the spin chain, and there is no other irreducible multiplet. 



4 The spectrum. 

Since all the irreducible multiplets are now identified, finding the spectrum consists in 
computing the action of the Hamiltonian on the highest weight vectors. The hamilto- 
nian (JI|) is su{2) invariant, hence it acts on the M-magnon subspace. In the magnon 
basis, this action is : 




The Hamiltonian act on the state of the form (||). Using eq. ( ^3] ) given in the Appendix, 
one realizes that the action induced on the polynomials ^{z) is : 

^ iV/ \p=l p<q Z Pq Z 1PJ 

= (Tr) EW-W) ( 20 ) 

^ ' p=l 

In Eq. (^0]) one recognizes the Calogero- Sutherland Hamiltonian at a special value of 
the coupling constant. In other words, the spin chain in the M-magnon sector has been 
mapped on the M-body Calogero problem. The last equality in (p0|), gives the energy 
of a multiplet {m p }: 

E({m p }) = Y: (^) 2 m p (m p -iV) 



5 Wave functions and statistics. 



Only the wave functions of the h.w. vectors are relevent since those of their descendents 
are obtained by recursive action of the transfer matrix. We now show how recent results 
on the Calogero models can be used to find explicit expressions for these wave functions. 
The latter are based on the construction of operators, which we denote by Am in the 
M-magnon sector, intertwining the Calogero Hamiltonian and the free Hamiltonian: 

M 

H A M = A M A with A = £ z p d Zp {z p d Zp - N) (21) 

P =i 

These intertwiners were defined in [[IT], [HJ. One of their definitions is the following 
Vandermond determinant of the operators D p : 

A M = £ sign(a) D^-_\.--Dl 2 D ai 

a perm. 

In this formula, it is understood that Am is acting on antisymmetric functions. For 
example, for two magnons: A 2 = Z\d Zx — z<id z ^ — + 



21-22 



The operators A^ are antisymmetric. Therefore, the symmetric wave functions are 
obtained by acting first with the antisymmetrizer on the plane waves z™ 1 ■ ■ ■ z M M , and 
then with Am'- 

*{z u ...,z M ) = A M (Det{z™% q ) (22) 

It is easy to check that the wave functions fl2"2p are symmetric polynomials vanishing at 
coincident points. If the rapidities are such that 1 < m p < (N — 1), these polynomials 
have degree less than (N — 1) and satisfy the condition (TH|). They thus are the wave 
functions of the h.w. vectors. In other words, since the plane waves z m are the wave 
functions of the one-magnon h.w. vectors, the operator T> = Am ° Det map tensor 
products of M one-magnon states into M-magnon states: 

\mi) ® • • • ® \m M ) — ► | {mi, ■ ■ • , m M }) (23) 

This map is a generalization of the Slater determinant in the sense that it implements 
the rapidity selection rules: if two of the rapidities m p and m q are either equal or differ 
by a unit, then the resulting wave function vanish identically. The fact that the result 
vanish if two of the rapidities coincide is obvious from the definition, while the fact that 
it vanishes if they differ by a unit results from an explicit check on the two-magnon 
case (which has all generality thanks to the symmetry of the wave functions and the 



recursive definition of Am given in 



6 The spinon description. 



The magnons are the excitations over the ferromagnetic vacuum; the excitations over 

thp flntifprrnmflcmptir vflpnnm arp rnnvpnipntlv Hpsrrihprl in tprms nf sninnns 



For N even, the antiferromagnetic vacuum corresponds to the alternating sequence 
of symbols 010101 • ■ -010. Its rapidities sequence is {mj = 2j — l}f=x. The excitations 
are obtained by flipping and moving the symbols and 1. We classify the sequence by 
their number M of rapidities. The spinon number N sp of a sequence is then defined 
by M = N ~^ sp . Since M is an integer, (N — N sp ) is always even. The spin S z of the 
Yangian highest weight vector of the sequence is S z = ~(iV — 2M) = ^p. 

A sequence of rapidities {mfj = 1, • • • ,M — N ~^ sp }, in the N sp sector, can be de- 
composed into (M + 1) elementary motifs. As we defined them in Sect.l, an elementary 
motif is a series of consecutive 0. We will think about the elementary motifs as the 
possible orbitals for spin half objects, which are called spinons. At fixed N sp , there are 
(1 + orbitals. To a spinon in the j th orbital, with j = 0, • • • , — we assign a 

momemtum k = ^jf. Thus, the spinon momemta vary from zero to ko = ^ ( N ~^ sp ^ . 

N-N av \ 



2 



By convention, a sequence of rapidities {rrij} corresponds to the filling of the (14 
orbitals with respective occupation numbers rij = rij + n~ = — rrij — 2), with 

m = —1 and tum+i = N+ 1. The length of the j th elementary motif is then Qj = rij + 1. 
By construction, the total occupation number is the spinon number : 



E (n+ + nj)=N sp (24) 

3=0 

Since an elementary motif of length Q corresponds to a spin representation of 
su(2), the full degeneracy of the sequences is then recovored by assuming that, at fixed 
spinon number, the spinon behaves as bosons. Notice that this property is specific to 
the su(2) spin chain. 

The spinons are not bosons but "semions" since the number of available orbitals 
varies with the total occupation number [TJ]. In particular, the spinons are always 
created by pairs. 

The energy of a collection of N sp = N — 2M spinons is : 

E-E vac = E (M) + ^^2(M-j)(M-^+jX ( 25 ) 
+ ££(M-sup(j,j"))nX; 

j,f era' 

with E (M) = \M(M - 1)(4M + 1) - (N - 1)M 2 . 

The low energy, low temperature, behavior is classified |15| as the level one su(2) 
WZW conformal field theory. In the spinon description the states consist of semi-infinite 
sequences of symbols and 1. The two primary states, which correspond to the two 
integrable representations of the su(2) current algebra at level one, are the vacuum, 
with sequence 010101 • • • and the spin half primary, with sequence 0010101 • • •. The 
excited states are given by finite rearrangement of the primary sequence. The Virasoro 
generator acts as L = J2j( m< j ~ m j), where {mj} are the primary sequences. Its spinon 
representation is : 

U=(^) 2 + T, .Unt + nj) (26) 



with N sp = J2j( n ~^ + n j), the total spinon number. The excitations over the vacuum 
possess an even number of spinons, wheras those over the spin half primary contain an 
odd number of spinons. 



7 The thermodynamics. 

Following ref. ||, the thermodynamics can be derived from the spinon description using 
methods similar to the thermodynamic Bethe ansatz [|T^] . 

First, we consider the system with a fixed spinon density D sp = In the iV — > oo 
limit, the spinon orbitals are labeled by momenta k continuously varying from zero to 
^° = W ( N ~^ sp ) = 7r(l — D sp ). Let n±(k) be the spinon occupation numbers of the k th 
orbital. By definition, they satisfy : 

E / ^n a (k) = D sp (27) 



In the continuum limit, the energy per unit of length is : 

E y-^ f k o dk 

— = £{D ap ; n±(k)) = S {D sp ) + V / — e (k)n a (k) 

iv ~ Jo An 



+ E/ ^V(k,k')n a (k)nAk r ) (28) 



-;, Jo (2vr) 



with £ (D sp ) = (|) 2 (1 - D sp ) 2 C-±f^), e (k) = §(fc - k)(k + k - tt) and V(fe, fc') = 
|(A;q — sup(/c, fc')). In each orbital the spinons behave as bosons, therefore the entropy 
of a configuration of n±(k) spinons is : 

S r k " dk 

- = 5(£) v ; n±(fc)) = £ / T~ (M k ) + X ) lo sK(*0 + 1) " M*0 logn„(A0) (29) 
iv " Jo 2vr 

The free energy per unit of length is : 

^- = T = E-TS-hM (30) 

with h the exterior mag netic field and M = J k ° §(n+(k) - n_(Jfe)) the magnetization. 

At fixed spinon density, the thermodynamic equilibrium is determined by minimizing 
(|50"D with respect to the variation of the spinon occupation numbers subject to the 
constraint (pTp. This gives : 

n ^ = e^l-A) _ x ( 31 ) 

where A is the Lagrange multiplier and e{k) is the dressed energy defined by : 

xf r k Q riy 

At fixed density, the Lagrange multiplier is determined from the constraint (p7|) . This 



The spin chain corresponds to a spinon gas of arbitrary density; i.e. the spinon 
chimical potential ll = -^j— is zero. Minimazing the free energy with respect to the 
density fixes A to be ft A = — log (2 cosh(/3/t)). The constraint (p7|) then gives the mean 
density D sp . 

The coupled eqs. (|3l|j3^) can be solved. Deriving twice eq. (|32|) with respect to k 
gives : 

d 2 e(k) 



dk 2 



(33) 



with the boundary conditions: e(k ) = e(0) = 0, and e'(0) = —e'(k ) = |. Eq . (|3~3"D is 
integrated by introducing the dressed momenta P = %■ It varies from — ~ to ~, and it 



satisfies ^ 



dk 2 



d(p 2 /2) 
de 



. As function of p, the occupation number are then given by : 



n± = exp [/3(?7 ± e dr ± /*)] 



(34) 



with 



Vip) = 2 



,7T, 



P 



and exp (/???) 



sinh(/3e dr 



sinh(/?/i) 

In the limit h — > 0, the dressed energy is ~ /ie^. 

The free energy is found by integrating the thermodynamical relations : 



(35) 



The result is the following simple answer : 



T = - T 



I dp 



71 



log 



smh(P(e dr + h)) 
sinh(/3/i) 



(36) 



(37) 



Notice that this is the free energy of a gas of non-interacting particles which, in the 
limit h — > 0, have energies given by r/(p). 



8 Appendix. 

In this Appendix we prove the eqs. ([10|JTT|) . First we compute the action of T (x) 

on these magnon states. We recall that T (x) can be written as T (x) = 1 + 

J2ij ( jzl) . -PT' w h ere is th e projector on the spin (a.,- = — ) acting on the j th spin 
only. The projector P~ on the M-magnon states (^j) gives states with one spin down 
marked of the form : 

Pr|¥) = |*i> = M £ Vw™*7^'--*»J n > (38) 

They corresponds to polynomials ^(z%\z2, • • • , Zm) symmetric in z 2 ,---,zm with the 
point z\ distinguished. To evaluate the action of J2ij{L n )ijPj~ , we remark that on this 
class of states, Ly acts as follows : 

V LMt) = \(LV)<) (39) 



with, 

M 

+ E 

(40) 

k q=l 

where 6jk = uji /{uji —u> h ), and K lq permutes the indices in position 1 and q. Therefore, 
\(L^)i) = Y!,j{L n )ijPr\^f) can be recursively computed. Then, 

E(^ B )^I*> = ^EI(^ B *)*> ( 41 ) 

is obtained by symmetrizing the wave function coefficients of |(L n \l/)j) in all its indices. 
I.e. its wave function coefficients, denoted (L n P^) nir .. nM are : 

(L n P*) ni ,., nM = |l + E K lp \ {L n ^) m .., nM (42) 

We now translate this action on the wave function coefficients into an action on the 
polynomials. We recall that in the basis of polynomials Qk in one variable of degree 
(N — 1) specified by Qk(w n ) = S^, the matrix elements of the derivatives are: 

iV + I u j 

{zd z -^-)Q j {z) = -E- - k Qk{z) (43) 

1 k+j u ~ u 

zd z {zd z -N)Q 3 {z) = (Q^-Qjiz)) 

The differentials D p , introduced in eq. ([12]) , acts on polynomials of the form @; i.e. 
they preserve the form of these polynomials. Moreover, by comparing the formula, we 
recognize in eq.(|40"|) the operator [l>i — ^p)- Since the operators D p are covariant 
by permutation of the coordinates, the polynomial of (L n P^) is given by acting with 
J2 P {p p — on ^ ■ Resuming all the contributions, we obtain : 

T-{x) = [ 1 + E : T A_n ) ^ 

\ p=l ' 2 J 

The action of T_ + (x) on these subclass of magnon states can be computed using the 
same methods. Its action is given by eq. (|TTp . It is remarkable that the differentials D p 
appear naturally in the study of the spin transfer matrix. 
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